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Abstract 

Jet broadening is an event-shape variable probing the transverse momenta of particles 
inside jets. It has been measured precisely in e + e~ annihilations and is used to extract 
the strong coupling constant. The factorization of the associated cross section at small 
values of the broadening is afflicted by a collinear anomaly. Based on an analysis of this 
anomaly, we present the hrst all-order expressions for jet-broadening distributions, which 
are free of large perturbative logarithms in the two-jet limit. Our formulae reproduce 
known results at next-to- leading logarithmic order but also extend to higher orders. 



1 Jet broadening 



Event shapes measure geometric properties of collider events. They are conceptually simpler 
than the more commonly used jet observables and have been measured with exquisite precision 
at lepton colliders such as LEP [U El El II]- By comparing to theoretical predictions, these 
results can be used to extract the strong coupling constant and to search for new particles. 
The interest in event shapes has been renewed recently after the fixed-order results at next- 
to-next-to-leading order (NNLO) in QCD perturbation theory have become available |6]. 
However, even with these corrections at hand, the perturbative uncertainty still dominates the 
error on the extracted value of a s . The convergence of the perturbative series can be improved 
by resumming those perturbative corrections which become dominant near kinematical end- 
points. This resummation was performed at next-to-leading logarithmic (NLL) accuracy for 
the event shapes measured at LEP; however, to improve the accuracy of a s it is necessary to 
also resum subleading logarithms of higher order. This has been achieved for the event-shape 
variables thrust [7] and heavy-jet mass [8]. In the present paper we derive, for the first time, 
all-order formulae for the total and wide jet broadenings, which are free of large logarithms 
and allow the resummation of higher logarithms also for these quantities. 

The jet broadening is defined as follows. With each hadronic event in an e + e~ collision 
one associates a thrust axis n^, which is the direction of maximum three-momentum flow. 
The particles of the event can then be divided into two groups: those moving in the forward 
hemisphere with respect to the thrust direction (p, -Ut > 0), and those moving in the opposite 
hemisphere. For simplicity, we assume that ut points to the left and refer to the two groups of 
particles as the left-moving and right-moving ones, respectively. The left-broadening is defined 
as the sum of the absolute values of the transverse momenta of the left-moving particles [91 [10] 

b L = \^\pt\ = \'52\PiX™r\, (!) 

and analogously one defines the right-broadening 5r. Usually, one normalizes the broadenings 
to the center-of-mass energy Q = y/s and defines Blr = bL,R,/Q, but we prefer to work with 
the dimensionful quantities bL,R- What is measured experimentally is the total broadening 
br = b L + bji as well as the wide broadening bw = max(6 i7 bu). 

In the following, we are interested in the region of small but perturbative broadenings, 
Aqcd <C &l ~ bR <C Q, where Aqcd is a typical scale associated with non-perturbative strong- 
interaction physics. Despite the fact that the broadenings are in the perturbative domain, 
fixed-order perturbation theory breaks down at small values of b^R, because large logarithms 
arise. For the fraction of events with total broadening smaller than br, for example, one 
obtains at next-to-leading order, 

R(B T ) = — [^dbT ^- = 1 + ^ (-4 In 2 B T - 6 In B T - 7 + vr 2 ) , (2) 
do J db T 2ir 

up to terms which vanish in the limit where Bt = br/Q goes to zero. Here and below we 
normalize to the tree- level total cross section cr - I n the result the Sudakov logarithms 
are manifest. To obtain reliable predictions for small broadening, these enhanced corrections 
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must be resummed to all orders. At leading double-logarithmic order this was achieved in 
To this accuracy the broadening can be written as a product of two jet functions in 
Laplace space. An improved version of this result, valid also at the single logarithmic level, 
was later presented in [12] . However, an all-order formula for broadening which is free of large 
logarithms was missing. Near the two-jet limit, event shapes such as thrust factorize into a 
convolution of a hard function, two jet functions, and a soft function, and this factorization 
forms the basis for an all-order resummation of logarithmically enhanced corrections. This 
was shown in [13] for a large class of event-shape variables, but it was also pointed out that 
this factorization breaks down for broadening. The same class of event shapes was recently 
reanalyzed in Soft-Collinear Effective Theory (SCET) in [13], which concluded that the usual 
effective-theory power counting breaks down for broadening. 

While the naive soft-collinear factorization indeed breaks down for broadening, we never- 
theless manage to derive in this paper all-order formulae for the total and wide broadening 
distributions, which are free of large logarithms. The reason is that the breaking of factoriza- 
tion has a very specific origin. In the effective theory, it manifests itself as a collinear anomaly, 
which generates an additional dependence on the large momentum transfer Q in the product of 
the jet and soft functions. In the effective theory the collinear anomaly is a quantum anomaly 
in the usual sense, that a symmetry of the classical (effective) Lagrangian is not preserved by 
the regularization. The factorization analysis is similar to the one for small-gx resummation 
in Drell-Yan production [15]. As in this case, in an intermediate step one needs to introduce 
additional regulators beyond dimensional regularization in order to obtain well-defined ex- 
pressions in the effective theory. When the regulators are removed in the final predictions for 
physical cross sections, the anomalous Q dependence remains. The regulator independence 
of the product of jet and soft functions gives a strong constraint on the dependence of the 
individual functions on Q, implying that this dependence must exponentiate [T5| IT6]. 

The analysis of jet broadening is complicated by the fact that not only collinear modes, 
but also a soft mode (using SCET n terminology), whose momentum components scale as 
p^ ~ bx, give a leading contribution to the cross section. In contrast to the present case, this 
soft mode does not contribute to qr resummation in Drell-Yan production at small transverse 
momentum, because the corresponding loop integrals are scaleless and can be omitted after 
proper regularization [T5]. This is no longer the case for jet broadening, since the radiation 
is restricted to one of the hemispheres. Interestingly, we verified that starting at two-loop 
order the ultra-soft momentum region p^ s ~ b\jQ also gives non- vanishing contributions to 
individual diagrams in the presence of the hemisphere constraint. However, as explained in 
|15j . in the sum of all graphs these contributions cancel as a result of the KLN theorem, and 
consequently the ultra-soft region does not contribute to the broadening. 

We have stressed that the jet and soft functions relevant for broadening are not well 
defined without additional regularization. Leaving this issue aside for the moment, the naive 
factorization theorem for small broadening has the form 

h^k-^fafaj^l*** (3) 

x J L (p L -b s L ,pi,/i) J R (b R -b 3 R ,p^,fj) S{b s L , b s R , -pi, - Pr -, h) . 
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Figure 1: A typical event with small broadening consists of energetic collinear partons in each 
hemisphere (blue lines) accompanied by soft radiation (red wiggly lines) . The total transverse 
momentum with respect to the thrust axis ut vanishes in each hemisphere. 



The convolutions over b s L R arise because the physical broadening is the sum of the collinear and 
soft broadenings. The definition of the thrust axis ensures that the total transverse momentum 
vanishes in each hemisphere, so if the left-moving collinear partons have transverse momentum 
Pl, the transverse momentum of the soft partons in the left hemisphere must be equal and 
opposite (see Figured]). 

The hard function H(Q 2 ,^) = \Cy{— Q 2 — ie,fi)\ 2 is the square of the quark vector form 
factor, and is known to three-loop accuracy [171 [18]. The quark jet function for the left-moving 
collinear partons is given by 

Mb, P ± ^) = ^ (2vr) d 5(fi ■ p x - Q) - p^) 

(4) 

x S { b - 2 E Wl) (°\x a (0)\X) (X\xfi{0)\0) , 

where n M = (1, Hr) is a light-like vector along the thrust axis, n M = (1, —itr) is its conjugate, 
and for simplicity we drop the subscript L on the variables of the jet function. The first 
two ^-distributions ensure that the produced jet X has the desired energy and that its total 
transverse momentum has a given value jr 1 . In the absence of soft-collinear interactions, the 
collinear SCET Lagrangian is equivalent to QCD and the collinear quark field x( x ) can be 
identified with x{ x ) = ^ W^(x) ip(x), where ip( x ) is the QCD quark field and W(x) a straight 
Wilson line along the n M direction from — oo to x (see e.g. [19] for more details). In our 
computation of the jet function, we will use the standard QCD Lagrangian and Feynman 
rules. The jet function J~r for the right-moving collinear partons is obtained by exchanging 
n M -H- fi^ in the above formula. The soft function is obtained as 

Xl,Xr 
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2 j ' / \ 2 

ieX L j€X R 

(5) 

Here S n and S n are soft Wilson lines extending along the and n M directions. The final 
states are split into left and right-moving particles, X = X^ + Xr, which contribute to the 
respective broadenings. 



3 



Instead of working with the cross section in momentum space, it is more convenient to 
discuss the Laplace-transformed cross section 

d2 ° = [°°db L e~ TLbL f°°db R e- TRbR — (6) 
dr L dr R J J db L db R 

It is furthermore beneficial to Fourier transform in the momenta pj; and p^, after which the 
factorization theorem fl3]) takes the form 

1 -7^- = ^f (d ' 2) H{Q\») [ d d ~>xi [ d d ~\i 

do dr L dr R J J (7) 

where J7l, 3 r , and S are the Laplace and Fourier-transformed jet and soft functions, e.g. 

/oo p Jd-2 ± 

db L e~^ J _i^ e HpW J L (b L ,pt,ri . (8) 

Since by rotational invariance the jet functions only depend on the moduli \xj; R \ of the 
transverse-position vectors, we can integrate the soft function over the solid angles associ- 
ated with these vectors. Introducing the dimensionless variables zl,r = 2\ x l, r \/ t l,r, the 
naive factorization formula then takes the final form 



1 d 2 a 



a dr L dr R 



POO POO 

H(Q 2 ,/j,) dz L dz R J L (T L ,z L ,lj) Jr(t r ,z r ,ijl)S(t l ,t r ,z l ,z r ,ijl) , (9) 
Jo Jo 



with 

S{tl,t r ,z l ,z r ,(i) = — ^ / dVt^_ 2 d^_ 2 S(r L ,r R ,xj;,x^,fx) , 

J J (10) 

J l>r {t, z,h)=N (2ir) d ~ 2 - ( y ) ^(r, x 1 , //) . 

The normalization factor 

A/" = " d ~ 2 = - (11) 

(2vr) d - 2 (4vr) 1 - e r(l-e) 1 J 

is chosen such that for a d = 4 — 2e dimensional space-time the soft function = 1 at lowest 
order. For the lowest-order jet functions, we obtain 

Jl#[t, z) - r2(1 _ e) (1 - z2) 3/2- e • ^ 
2 Evaluation of the jet and soft functions 

As mentioned above, the Feynman diagrams contributing to the jet and soft functions are 
not well defined individually in dimensional regularization. To evaluate them, we introduce 
additional analytic regulators [201 EI] in the QCD diagrams contributing to the broadening. 
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Figure 2: Next-to-leading corrections to the jet function. The one-loop virtual corrections are 
scaleless and vanish. 



This is done as follows. Let us assume that the QCD vector current produces a quark which 
ends up in the left-moving jet, while the anti-quark ends up in the right-moving jet. We now 
raise all propagators along the quark line to a fractional power, 

_L_ _> K 2 )" (13) 

and analogously for the anti-quark line, but with a regulator (3 and an associated scale v^- 
Since the QCD diagrams are well defined without additional regulators, this operation is trivial 
for the full theory, but with the analytic regulators in place the effective-theory diagrams are 
now also well defined. The individual jet and soft functions in the naive factorization formula 
fl9]) exhibit poles in the regulators a and (3, which only cancel in the product of these functions 
entering the formula for the cross section. As is explained in detail in [JS], the regulator 
on the anti-quark line regularizes the Wilson line in the collinear jet function of left-moving 
particles. The usual eikonal Feynman rule for the emission of a left-moving collinear gluon 
with momentum k from the Wilson line in the SCET quark field \ is modified to 

(ylf n»n-p R 



n-k {n-kn-p R ) l+f} 

where pr is the total right-moving momentum, which fulfils n-pR = Q at leading power in br/Q- 
Similarly, the rule for the emission of a right-moving collinear gluon from the corresponding 
Wilson line is obtained by replacing z/ 2 — > ui, (3 — > a, n^ <->■ n M , and pr — > pi in the above 
result. Note that the usual eikonal identities no longer hold in the presence of the analytic 
regulators. The Feynman rules for multiple emissions from the regularized object are thus 
more complicated than for standard Wilson lines, but they can easily be derived by considering 
the corresponding regularized QCD diagrams. 

The regulators a and (3 also regularize the soft emission diagrams in the effective theory. 
The eikonal factor relevant for gluon emission from the regularized soft Wilson line S n is 
replaced by 

n M (vf) a n ■ p L 



n-k (n-kn- 

The same expression with V\ — >■ u 2 , a — > n^ 1 <H- and pl —> Pr holds for Sn, the Wilson 
line describing soft emissions from the anti-quark. 

With the additional regulators in place, we may now evaluate the next-to-leading order 
corrections to the jet and soft functions. The Feynman diagrams relevant for the jet function 
are shown in Figure [2J Let us first consider the simple case pj; R = 0, for which the one- loop 
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expression for the jet function can be obtained analytically. In terms of the MS coupling 
constant a s = a s (fi), the next-to- leading order result for the bare diagrams relevant for the 
left-moving jet reads 



J L (b, P ± = 0) = S(b) 



C F a s 



2tt rn 



^lE 



1 b \b) 



2c 



X 



(1 



4r(2 + a) r(q - g) (v\_ 
T(2 + 2a-/3) 



^2 



(16) 



The first diagram in Figure [2] produces the first term inside the brackets, which is well defined 
without analytic regulators and thus has been evaluated setting a = (3 = 0. The second and 
third diagrams give identical contributions and produce the second term inside the brackets. 
These graphs are only well defined if the analytic regulators fulfill a ^ (3. One can thus only set 
one of the regulators to zero. Expanding then in the second regulator one encounters a pole. 
The fourth diagram vanishes. Note that at next-to-leading order the regularized jet function 
is gauge invariant, since only real-emission diagrams contribute. At higher orders, the two jet 
functions and the soft function will likely not be gauge invariant individually, since analytic 
regularization breaks gauge invariance; however, gauge invariance will be recovered in the 
product of the three functions. The result for the jet function J7r follows from the expression 
for J L given above by exchanging a •<-> /3 and i>\ •<->■ i>2- For p^ p> 7^ 0, the computation of 
the jet functions becomes rather non-trivial. Restricting ourselves to the divergence in the 
analytic regulators, we find up to terms of (9(e) 



J L (b,p J 



8(b- P - 



2C F a s 



71 



a 



P b \bJ 



2c 



462 



-l-2e 



Q 7 



(17) 



where p = \p L \ , and the dots represent non-singular terms in the analytic regulators. Perform- 
ing the Fourier and Laplace transformations, the bare jet function takes the form 



J l (t,z) = J { °\t,z) 



X 



C F a s 



71 (3 



i,(i +1 » ( AV 2 ^)(*T(|) , + 



5) 

where we have defined f = re 7B , and 3^ (r, z) is the tree-level jet function from f|T2|) . 

The soft function is obtained by evaluating the diagrams in Figure [31 The corresponding 
analytically-regularized Wilson-line integrals are simpler to evaluate than those for the jet 
function, so we are able to obtain a closed expression at one-loop order. In momentum space, 
we find (with pi = \Pl\) 



)5 d -\pi)5 d -\pi) + ^{^y 



7T 



2 \ a / , ,2 \ P 



QJ \Q 



f3 — a 



p- L M 5(b L - 8(b R ) 5 d ~\pi) -{L^R) 



(19) 
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Figure 3: Next-to-leading corrections to the soft function. The one-loop virtual diagrams are 
scaleless and vanish. 



The Laplace and Fourier transforms can be performed explicitly, yielding 



S(t l ,t r ,z l ,z r ) = 1 + 



C F a s 1 2 1 ""-^- 2£ r(-q -P-2e) 
ii f3 — a 



ni 



w ( | 



2 \ a / , ,2 \ 



a+/3+2e p 



a + /? + 2e l-a-/3-2 £ 1 2 , 
S > n > 1 - e, - (£ #) 



(20) 

We now expand this expression by taking the limits /3 — >■ 0, a — > 0, and e — > 0. The order 
in which the analytic regulators are taken to zero is arbitrary, but it is important that the 
limit e — >■ is taken at the end, since only then the QCD result is independent of the analytic 
regularization. The final expression for the bare soft function is 



S(t l ,t r ,z l ,z r ) = 1 + 



4tt 



a Q 



2 In 



(21) 



+ 8 Lis 



1 + zl 



zl + 1 



4 In' 



zl + 1 



5tt 2 



(L R) 



where tl = TLe lE . The coefficients of the 1/a poles are equal and opposite to those in the jet 



functions, see (TT8|) . so that these divergences cancel in the product JlJrS. 

Let us now use the expressions derived above to compute the differential cross section at 
one-loop order. In this approximation we only need the convolutions of the tree-level soft 
function with the one-loop jet functions and vice versa. Since the tree-level soft function 
involves 5-functions in the transverse momenta, we only need the jet function J^L{b,p ± = 0) 
given in ([TBI) . The resulting expression can be refactorized in the form 



d 2 a 



cr db L db R 



H(Q 2 )Z(b L )Z(b R ) 



where 



£(&) = 5(b) + 



C F a, 



4tt Til 



2c 



6-2e 



(22) 



(23) 
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The dependence on the analytic regulator has indeed canceled among the jet and soft functions, 
and the anomalous logarithmic dependence on Q is now manifest. The above expression for 
S(6) is valid for arbitrary values of e. Expanding around e = 0, and carefully treating the 
resulting distributions in b, one finds the one-loop divergence 

=(»> = ^T «(»>(!-! I"? + f)+0(0- (24) 

This is equal to minus one half of the divergence of the bare hard function, so that the full 
cross section (122]) is indeed finite. In terms of renormalized objects, this result implies that 
the convolution of the jet and soft functions has the same anomalous dimension as the hard 
function, such that the cross section is renormalization-group (RG) invariant. We have used 
expression (j22p to compute the leading singular terms for small bx in the total broadening and 
reproduce the one- loop result (T5J). Having performed these one-loop checks, we now turn to 
the resummation of the logarithmically-enhanced corrections. 



3 Resummation 

The explicit results in the previous section show that the jet and soft functions contain diver- 
gences in the analytic regulators. These divergences cancel in the product of these functions, 
but they leave behind large logarithms of the momentum transfer over the broadening, which 
should be resummed to all orders in order to get reliable predictions. 

Interestingly, the one-loop divergences in the analytic regulators are multiplicative, such 
that the product 

P(Q 2 , T L , T R , Z L , Zr, fl) = J L (t L , Z L , fi) Jr(t r , Zr, fl) S(t l , t r , z l , z r , fl) (25) 

is finite even before the integrations over zl and zr in (j^J) are performed. To show that this 
property holds to all orders in perturbation theory, we now study a modified event shape, in 
which we consider the emission of two soft photons with momenta p lL and p lR in addition to 
the QCD partons. We then determine the thrust axis using all final-state particles, including 
the photons, but measure the broadening of the hadrons only. In the presence of photons in 
the hemisphere, the total hadronic transverse momentum no longer adds up to zero, but will 
be equal and opposite to the photon momentum. 

The emission of soft photons is described by QED Wilson lines, analogously to the ones 
appearing in the QCD soft function in ([S]). For our purposes, it is sufficient to treat the 
photons classically and to consider the case where each hemisphere contains one photon with 
fixed momentum. The QED soft function is then given by 

Sqed{P~, l ,P~, r ) = 4e = = , (26) 

n-p lL n-p lL n- p lR n ■ p lR 

and it is just an overall factor multiplying the hadronic cross section. We now analyze the 
differential cross section 

1 d 8 a 



oc #(Q 2 ,/i) Jdb s L J db s R J d d ~ 2 pi jd d ~ 2 pi 



^qed(p 7l ,P 7J j) db L db R d 3 p lL d 3 p^ R J J J J (27) 

x J L (b L - b s L , -pi, n) J R (b R - b s R , -pk, n) S(b s L , b s R ,pi - P^ l ,Pr- - p^ R , fi) , 
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which fulfills a factorization theorem involving exactly the same jet and soft functions as in 
()3]), except that the soft and collinear transverse momenta now add up to minus the photon 
transverse momenta instead of zero. Taking the Laplace transforms in the broadenings, and 
Fourier transforming with respect to the photon transverse momenta p^ L and p^ R , we find 
that the cross section is given by the same product of the Laplace and Fourier-transformed jet 
and soft functions as in (j7j), but without the integrations over xj; R . Since (127)) gives a physical 
cross section, it follows that the product P of the renormalized jet and soft functions in (125)) 
must be well defined and independent of the analytic regulators to all orders in perturbation 
theory. On the left-hand side of that equation we have indicated that this product carries an 
anomalous dependence on the hard scale Q. The naive factorization formulae ([3]), ©, and 
(1271) thus do not achieve a proper scale separation. 

The fact that the product P is independent of the regulator scales v\ and v 2 has important 
consequences, since the two scales enter the various functions in a very specific way. To derive 
these, let us send /3 — > first. As long as a is non-zero the corresponding limit is regular 
and the jet and soft functions are z/ 2 independent. What remains in the individual functions 
is a dependence on the regulator scale v = v\. However, the v independence of the product 
implies that d In P/d In v 2 = 0, where 

2 2 — 

InP = \nJ L (\nv 2 fl; r L ,z L ,fj) + In J R (ln T R ,z R ,fij + lnS^ln^p 1 ; t l ,t r , z l , z r , fij . 

(28) 

We indicate the explicit form of the v dependence of the three functions in the first argu- 
ment. The use of tl rather than f R in the first argument of the soft function is a matter of 
choice. In the collinear sector of left-moving particles, the logarithmic dependence on v can 
for dimensional reasons only appear in the form of ln(z/ 2 f|), since no other momentum scale 
is introduced by the a regulator in (TT3|) . In the collinear sector of right-moving particles, on 
the other hand, the induced dependence on v is accompanied by n ■ p L , because the collinear 
Wilson line is regularized in analogy with ( I14p . By Lorentz invariance, a logarithm of the 
"foreign" momentum component n ■ pi can only appear in the form of n ■ p R n ■ pi = Q 2 , so 
that we encounter ln(V 2 /Q 2 ). Finally, in the soft sector the relevant scale is always given by 
the geometric average of the large and small scales, yielding \n(u 2 fL/Q) in the present case. 

We now call }'l{.L; Tl,z l ,{i) = \nJ L {L\ t^,z l ,[i) and similarly for the right jet function, 
and denote by fj/{L\ r L ,z L ,n) the n-th derivative of this functions with respect to its first 
argument. The condition that the sum ( 12 8 p as well as its derivatives with respect to \wQ must 
be independent of the analytic regulator scale v yields the conditions 

( In v 2 f I- r L , z L , /i) = (2- 1 - 1) f R n) (in ^; r R , z R , //) ; n > 2 . (29) 

It follows that the second derivatives of the two functions must be equal to the same constant 
k 2 , and all higher derivatives vanish: 

f L (L; r L , z L , n) = L 2 + k[ (r L} z L: /i)L+ k^(r L , z u \i) , 

(30) 

k (y(i) 

f R (L; t r , z r , (j) = L 2 + fcf (tr, z r , fi) L + k^(r R , z R , fx) . 
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We can now evaluate (1281) with the particular choice v 2 = Q/tl, for which 



InP = k 2 (jj) ln 2 (Qr L ) + kf (r L , z L , /i) - k^(r R , z R , ji) \yi(Qt l 



(31) 



where the dots represent Q-independent terms. Finally, the fact that the result must be 
left-right symmetric implies that k^(r,z,fi) + k 2 (fi) ln(/ir) = — k^(r,z,fi) — k 2 (fi) ln(/rf) = 
—2Fb{t, z, fi), and hence the final answer can be written in the form 



InP 



^\n 2 {Q 2 f L f R ) -F b (t l ,z l ,ii) \^{Q 2 ff)-F B {r R ,z R ^) \n{Q 2 f 2 R ) 

{62) 

+ lnW{r L ,T R ,z L ,z R ,fj,) , 



where the remainder function W is independent of Q and left-right symmetric. 

We can gain further information by exploiting the fact that the cross section (1271) and hence 
the product H(Q 2 , fi) P(Q 2 ,tl,t r , Zl, z r , /i) must be RG invariant. From the RG equation 
for the hard function [To] 



d 



dln/j, 
it then follows that 

d 



O 2 

2T cusp (a s ) ln\+4 7 9 (a, 
fi 2 



H(Q 2 ,v) 



(33) 



dlnfi 



Mm) = o, 



d 



dlnfi 



F B (r,z,fi) = r cusp (a s ) 



d 



(34) 



dln/j, 



W{t l ,t r ,z l ,z r ,h) 



2r cus p (o? s 



) In (h 2 t l t r ) - 47 9 (a s ) W(t l , t r , z l , z r , fi) . 



The first equation implies that k 2 must be a constant. From the fact that at tree level this 
constant vanishes it follows that k 2 = to all orders, since there would be no way to compensate 
the scale dependence of the coupling a s (/i). Next, using our explicit one-loop results in (TTHj) 
and (ED, we find that 



F B (r,z,fi) 



C F a s 



71 



ln(/!T) + In 



VT+z^+i 



0(a 2 s ) . 



(35) 



We are now in the position to state the main result of this paper, which is the corrected, 
all-order generalization of the naive factorization theorem (Q for the Laplace-transformed 
double-differential cross section. Relation (I32p . combined with the fact that k 2 = 0, implies 
that the anomalous dependence of the jet and soft functions on Q exponentiates, and that the 
cross section can be refactorized in the form 



d 2 a 



cr dr L dr R 



poo poo 

H(Q 2 ^) dz L dz R (Q 2 f 2 L )- FB{TL ^ (Q 2 f 2 R ) 
Jo Jo 



(36) 



x W(t l ,t r ,z l ,z r ,ii) 



This result contains two sources of Q dependence: one arising from the hard function H(Q 2 , /j,), 
and an additional one stemming from the collinear anomaly. Once the RG equation ( 1331) for 
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the hard function H(Q 2 ,fi) is solved and the result is evaluated at a scale \i ~ 1/tl ~ l/ r #; 
the above formula no longer contains any large logarithms in the perturbative series for the 
functions Fr and W . 

It is convenient to factor out the leading-order jet functions from the remainder function 
W by rewriting 

W(t l , t r , z l , z r , /i) = — ' . - Zb ' , W(t l , t r , z l , z r , fi) , (37) 

(1 + z£) (1 + 4) 

where IV = 1 + 0(a s ). At NLL order, we can evaluate the cross section using the tree-level 
result W = 1 and the one- loop expression for F B given in ( 131)]) . This yields 

H(Q 2 , /i) Outx)-^ O'ts)-'- Jfo) I( m ) , (38) 



cr dr L dr R 



where 



and 



VL,n = In (Q*fl R ) , (39) 

7T 



/■°° Z /,/1 + z 2_l 4r? 

m = L dz jT^r^r^) =r^^i+".2+,,-i). (40) 

The function /(r/) was called (2/X) v in [12]. Finally, the relevant expression for the hard 
function at this order reads [22] 



Pq 



' 1 - - - lnr ) + ( K - ^ ) (1 -r + lnr) 



V r / V Po, 

(41) 

Pi i 2 , 3/?o , 
H — In r H in r , 

2f3 2 J' 

where r = a s (^)/a s (Q), (3 and /3i are the first two expansion coefficients of the QCD f3- 
function, and K = (y — ^) Ca — y TpUf. At NLL order we can further approximate 

_ C F a s (n) Q 2 
Vl,r ~ t? = In ~T ' ( 42 ) 

since hi (/rf^) is a small logarithm and counts as 0(1). Then the cross section exhibits only 
simple power dependence on tl^r, and the Mellin inversion from Laplace to b^ R space can be 
performed analytically. We finally obtain 

L^W^U^'U!*)'^). ,43) 
cr db L db R T 2 (r]) b L \ fi J b R \ fi J 

For the total broadening br = bi + bR and the wide broadening bw = max(&£, &r), we find 



H (Q >/*) T^rr t — 1 (v) 



o"o dbx ' r(2r/) b? \ [i 

1 rfcr 2 2r/e- 27Bf ' 1 /fr 



2'/ 



(44) 



cr ao^y T 2 (l + rj) b w \ \i 
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Figure 4: The upper panels show the LO and NLO coefficients in the expansion (l4"6"j) of the total 
broadening cross section as a function of L = In ( bj< /Q), as computed by the generator EVENT2 
[23] . The lower panels show the difference of the full results to the NLL approximations. 



In these results the resummation of the LL and NLL terms is accomplished by setting /i ~ 
~ b R or fi ~ bx,bw, respectively 

Our expressions for the cross sections in ( T44")) are completely equivalent to the NLL results 
obtained in [12], while an earlier result derived in [11] misses the factor I 2 (r]) and is thus only 
correct at leading double-logarithmic order. In order to relate our expressions to the ones 
derived in these papers, we choose fi = b T and rewrite the coupling constant a s (br) at the low 
scale in terms of a s {Q) using 



a s (b T ) = TT - 



A ln(l + ") a s{Q) 
/3 (1 + w) 4tt 



u= ^m ln bi (45) 

47T Q 2 



As a further test of our calculations, we compare our results to the fixed-order calculation of 
the jet broadening cross section. To our knowledge, an analytical result for this quantity is not 
available in the literature - not even at 0(a s ), which is leading order (LO) for the distribution. 
Numerically, the cross section can be computed to next-to-leading order (NLO) using the 
EVENT2 generator [23J. The NNLO result has been obtained in [51 IS]- Figure H] shows the 
difference between the NLO result as computed using EVENT2 and the NLL singular terms 
predicted by our resummation formula. In the figure we show 

bTda a ^A(b T ) + (^) 2 B(b T ). (46) 



cr db T 2n \ 2ir 
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The leading singular terms for small bp, which are predicted at NLL and can be obtained by 
expanding (|44j) . are given by 



A NLL (b T ) = C F (-8L - 6) 



B NLL (b 



'T 



IL 3 + 72L 2 + ( 92 - - 641n 2 2 



(47) 



+ C f Ca 



" 88 T 2 


/4tt 2 




— L 2 + 






. 3 


V 3 





CpTpUf 



3 



- L 
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with L = ln(6^/Q). At leading power in bp/Q, all terms in A{bp) are predicted by our results, 
while the constant term in Bibp) can only be obtained by extending our analysis to NNLL 
order. For b? — > 0, the singular terms dominate the cross section. We should thus find that 
the difference AA = A{b T ) - A NLL (b T ) tends to zero for b T 0, while AB = B(b T )- B NLL (b T ) 
approaches a constant in the limit of small broadening. This is indeed what is observed in 
Figure HI Let us note that a similar numerical comparison was also shown in [12] . However, 
we have run EVENT2 with a very low cutoff of 10~ 16 and have computed 10 8 events using 
quadruple precision, which was computationally not yet feasible at the time when that paper 
was published. 



4 Summary and outlook 

Using effective field-theory methods, we have presented the first all-order factorization theorem 
for jet broadening distributions in e + e~ collisions, in which all large logarithmic corrections 
arising for small broadening are resummed. Our result has several interesting features. First, 
the broadening receives both collinear and soft contributions, with the collinear jets in each 
hemisphere recoiling against the soft radiation. These recoil effects were discovered in [12] 
and become non-trivial at NLL order. Second, the individual jet and soft functions are only 
defined if additional regulators beyond dimensional regularization are introduced. We have 
used analytic regulators for that purpose. The individual functions then contain divergences in 
these regulators, which cancel in the product of jet and soft functions entering the formula for 
the differential cross section. When the regulators are removed an anomalous Q dependence 
arises, which constitutes a new source of large logarithms. Such an anomaly also appears 
in transverse-momentum resummation in Drell-Yan production [T3] and in electroweak Su- 
dakov resummation [16]. Using similar arguments as in these papers, we have shown that the 
associated anomalous Q dependence of the broadening distribution exponentiates. 

The analytical results obtained in this paper allow us to resum logarithms to NLL accuracy, 
and in this way we have reproduced the results of [12] • In addition, we have checked the 
resummed expressions against the numerical fixed-order results and find that our expressions 
correctly reproduce the numerical results for small broadening. It would be interesting to 
obtain NNLL accuracy for the total and wide broadening and to combine the resummed 
expressions with the existing fixed-order results at NNLO obtained in [51 16]. This could be 
used for a precise determination of a s from the existing measurements of the total and wide 
broadening distributions. For the event-shape variable thrust, the resummation has been 
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performed to N 3 LL accuracy [7], and this result was used for a precision determination of a s 
in [21]. The resulting value is much lower than the world average of a s , and it is important 
to cross check whether other event shapes lead to a consistent value, if only to gauge whether 
the uncertainty estimates are reliable. To obtain NNLL accuracy for broadening, one needs 
the one-loop expressions for the jet and soft functions and the two-loop expression for the 
anomaly function Fb- The one-loop soft function has already been given in this paper. The 
computation of the one-loop jet function is however more complicated. Finally, to obtain 
the two-loop anomaly coefficient one needs to compute the two-loop divergence of the soft 
function in the analytic regulators. We hope to report results for these quantities in a future 
publication. 

Note added: While this paper was in preparation the preprint [2S] appeared, which also 
considers the resummation of jet broadening in the context of SCET. This work contains 
the factorization theorem ([3]), but uses an alternative analytic regulator for the jet and soft 
functions. In contrast to our work, the regularization is introduced directly in SCET instead 
of QCD, and it is not clear to us whether QCD is recovered once this regulator is removed. 
More importantly, the result for the total broadening presented in [25] does not account for 
the recoil effects discussed in [T2] and is therefore missing the factor I 2 (r]) in (1431) and (JHJ). It 
is thus incorrect beyond the leading double-logarithmic order. 
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